DYNAMICS OF SINGLE CHAINS OF SUSPENDED FERROFLUID PARTICLES. 

S. Cutillas and J. Liu, California State University Long Beach, Department of Physics and Astronomy, 

1 250 Bellflower Boulevard, 90840 CA, Long Beach. 


Introduction 


Super paramagnetic particles is a new class of 
magnetic materials that exhibit a susceptibility X near 
the value of one while most of natural paramagnetic 
materials have a magnetic susceptibility at least five 
orders of magnitude lower. Ferrofluid emulsion is one 
of such systems. A collection of iron oxide nanoparti- 
cles (in our case Fe 3 0 4 with a radius of 6nm) dispersed 
in kerosene is called a ferrofluid which has various 
properties [11. The most important is that each nano- 
particle has a permanent dipole moment. Due to the 
small size of the nanoparticles, they only have a single 
domain magnetization. The thermal agitation random- 
izes the iron oxide particles so that the total magneti- 
zation of the ferrofluid is equal to zero when no mag- 
netic field is applied. Once an external magnetic field 
is applied, the nanoparticles are oriented slightly in the 
field direction and the ferrofluid can have a strong 
magnetization. Particles, having a submicronic size, (in 
our case 0.23pm as radius) can be made of this ferro- 
fluid. If these ferrofluid droplets arc dispersed in wa- 
ter, some surfactant is added to prevent aggregation, 
the system thus produced is called a ferrofluid emul- 
sion. 

When ferrofluid particles are suspended in water, 
they diffuse under Brownian forces. Applying a mag- 
netic field to such a suspension leads to dipole-dipole 
interactions between particles. The amplitude of this 
interaction is described by a dimensionless parameter 
X, which is the ratio of the magnetic dipolar energy 
over the thermal energy [2]: 

18 kT 

Where a is the particle radius, // q is the vacuum 
permeability, H 0 the applied magnetic field, k the 
Boltzmann constant and T the absolute temperature. 
When X is greater than one, magnetic interactions 
dominate. Because of the strong anisotropy of the di- 
polar interaction, particles aggregate to form chains in 
the field direction (cf. Fig.l). 
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Figure 1: Magnetic field effect on a ferrofluid emul- 
sion. 


When the volume fraction is high, bigger aggre- 
gates are formed by the lateral coalescence of several 
chains [3,4,5]. While many studies and applications 
deal with the rheological properties [6] or the field- 
induced structural changes [3-7], this work is to better 
understand the chain dynamics. We will show, in the 
first section, that the kinetics of chain formation is a 
diffusion-limited aggregation. The second and the 
main part of this work will talk about measurements of 
the effective diffusion coefficient by dynamic light 
scattering (DLS) [8]. In particular, our experiments arc 
able to measure two different motions. One is from the 
center of mass of the chain and the other is from inter- 
nal motions (particle fluctuations). This is an extension 
of our earlier work [9,18]. 

Experimental conditions 

Our sample consists of ferrofluid particles dis- 
persed in water. The particle volume fraction is 10 5 to 
avoid multiply scattering as well as chain-chain lateral 
aggregation. The emulsion is introduced in a test tube 
of 2cm of diameter with care to avoid dust pollution. 



Figure 2: Experimental setup of dynamic light scatter- 
ing. The onset specifics our scattering geometry, with 
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chains perpendicular to the scattering plane, and the 
definition of q is also indicated. 

The tube is placed in the sample holder with an in- 
dex match bath to reduce reflections. The sample 
holder is then put inside a pair of Helmholtz coils as 
shown in Fig. 2. 

An Ion Argon laser beam (Ao=514.5nm) is focused 
through a lens (with a focal length equal to 25cm) to 
the sample. The detector is placed at an angle 9 from 
the initial path that specifies the momentum transfer q , 
defined in the experimental setup shown in Fig. 2. The 
scattered light is collected by a photomultiplier tube 
and analyzed by a digital correlator to give the inten- 
sity autocorrelation function. In self beating mode, 
often called homodyne mode, the intensity autocorre- 
lation function that the correlator measures 
<I(q,t)I(q,0)> is related to the dynamic structure factor 
S(q,t) with the relation: 

< f (q ,t) I (q ,0) > — < I 2 >= ccS 2 (q ,1) (2) 
And 

S(<7>0 = 2Xexp{i<7(r„ (t)-r p (0))} (3) 

V n ,p 

</ 2 > is the mean square amplitude of the intensity 
and a is a parameter that only depends on geometrical 
factors. V is the scattering volume and r n (t) is the posi- 
tion of the particle number n at time t. In general, in- 
terpretation of DLS experiments requires a model for 
S(q,t). 

In the absence of magnetic field, the suspension is 
isotropic without any interactions between particles 
and the structure factor reduces to a single exponential 
S(q,t)oc exp(-q 2 Dot). D 0 =(kT. ’/£) is the diffusion coef- 
ficient of a particle. j] a is the Stokes frictional 

coefficient of a spherical particle in a fluid of viscosity 
rj . If interactions or polydispersity can not be ig- 
nored, an effective diffusion coefficient has to be in- 
troduced. Formally, D^is defined as: 

D eff =— “T ln(5(<7,,)) 

q vt ,_> 0 

It reduces to D 0 if no interactions between particles 
arc present. Therefore, without magnetic field, we can 
use DLS to measure particle size. The particle radius 
was found to be a— 0.23/1 m with 7% of polydispersity 
through a fitting procedure called cumulant expansion 
[8]. In this case, if we vary the scattering angle 6 
(and so q) we do not have any change in the measured 
diffusion coefficient: it is ^-independent. 

When a magnetic field is applied, the system con- 
sists of particles and chains. Equation 2 can still be 
used to fit our data if we assume S(q,t)~exp(-C2 t) and 


Q =q 2 D eff . Q measures the frequency of scatters 
(particles, chains). 

The magnetic field is controlled by a computer 
(with a digital-analog converter board) connected to a 
remote power supply. The field was first calibrated 
without any magnetic material and a linear relation 
between the current and the field was found. Each 
Amp produces 67.0±0.5 Gauss. The relation giving X, 
for a given magnetic field strength is not simple since 
the susceptibility varies with the magnetic field. How- 
ever, an order of magnitude can be found with the 
equation X=0.01B 2 , where B is the external magnetic 
field applied, expressed in Gauss. For 200Gauss, 
X=400. Our expressions of X, given in this paper, are 
rigorous solutions of the equations of the magnetic 
field in heterogeneous media. 

Results 

I) Kinetics of chain formation 

When a strong magnetic field is applied (X=406), at 
a fixed q , we first notice a monotonous decrease of the 
effective diffusion coefficient with time. Physically, 
this means that chains arc becoming longer and longer 
with time (see Fig.3). 



Figure 3: Normalized effective diffusion coefficient 
versus time. (X=406, q= 4.2x1 0 6 m ! ) 

To obtain chain length from D eff , a theoretical 
model of rigid chains of N spherical particles is used 
[11,12] since it agrees very well with a recent experi- 
ment [10]. The result is shown below: 

D chain =f(N)D 0 

fff.-l 3I "< 2I ''> + 1 - 254 (4) 

J K } AN 

Here, we assume that the effective diffusion coeffi- 
cient measures the average motion of isolated chains. 
In very dilute system like the one we used, this is cer- 
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tainly true. From equation 4, we are able to get the 
mean number of particles per chain, N } from the effec- 
tive diffusion coefficient. 

This is illustrated in the Fig.4 where the effective 
diffusion coefficient is normalized by that of a single 
particle. Figure 4 shows a Log-Log plot of the results 
based on measurements of Fig. 3. 

Calculating the value of A as a function of the time 
gives us information about the kinetics of aggregation 
(i.c. /V versus time). 



Figure 4: Mean number of particle per chain versus 
time, kinetics of chain formation. 

Fitting the data with a power-law relation gives an 
exponent of 0.47±0.05. Thus, N depends on time as a 
square root (within our experimental uncertainties for 
this short chain formation). As expected, since the vol- 
ume fraction is very low, the kinetics of particle ag- 
gregation into chains is a diffusion-limited aggrega- 
tion. This is consistent with our earlier data and other 
work [9,18,131. 

II) Dynamics of single chains 

If we apply a magnetic field of 200 G to the sample 
for 6 hours, we can get chains long enough (more than 
19 particles) and far away (more that 70 particle di- 
ameters) that the chain growth slow down. Now, we 
can do other experiments. In this case, the absolute 
time is no longer important since chain length is al- 
most constant during the next 30 minutes. First, we 
vary 6 (or q) at a constant magnetic field strength, 
and measure the effective diffusion coefficient. 

Varying q allows us to select a window within 
which different length scales of motion will be probed. 
This can be seen in another way: q defines a charac- 
teristic length (1=2 7t /q) over which dominantly mo- 
tion will contribute to the intensity correlation func- 
tion. To be able to probe different length scales by 


varying the scattering angle is an important property of 
DLS. 

In this work, the experimental range of the scatter- 
ing angle is 5 0 < 0 < 130°. The equivalent range for 
the scattering wave vector is 1.4xl0 6 <<y(m 1 )<2.9xl0 . 
Thus, the characteristic length scale will have the 
range 0.2 1< / Qim) <4.42. In terms of the particle ra- 
dius: 

0.9 a < l < 20a . (5) 

We can see here that motions will be probed over 
quite a broad range in a single experiment from one 
particle radius up to twenty panicle radii. 

Next, we vary the magnetic field intensity (or X) 
and repeat the q-dcpendence experiment. 

Figure 5 shows four series of experiments per- 
formed at different magnetic field or coupling constant 
X (indicated in the onset). To work with the same chain 
length, these four experiments are done consecutively 
after the 6 hours at ^=406. The first one is for ^=406, 
the second one is for X^=17, the third one is for ^=157 
and the last one is for X=46. Each experiment takes 
around 30 minutes. This procedure allows us first to 
decrease the magnetic energy of interaction between 
chains to slow down again the kinetics of aggregation 
and second to cool the coils which were warm after a 
long-time experiment. At each fixed l value, we 
measure the effective diffusion coefficient D eff as a 
function of q. D eff is normalized by the single particle 
diffusion coefficient D 0 , and l=2n fq is normalized by 
the particle radius a. 



Figure 5: Dependence of the normalized effective dif- 
fusion coefficient versus the normalized scattering 
wave vector for different X values. 

Figure 5 shows a few interesting features. First, for 
all X values, D^is linearly proportional to qa. Second, 
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as X increases the slope of D eff versus q decreases. But 
all these straight lines converge to the same point 
(D eJ /Do= 0.16) as q approaches zero. 

To understand the q dependency of D eff , let us look 
at the probing length used. In our study, the character- 
istic length scale probed is l=2n /q which is in the 
range of 0.9<Z/ci<20. When l is much larger than the 
particle radius ( l/a -» °°) we are mainly sensitive to the 
center of mass diffusion of the chain. D eff is then the 
diffusion of the entire chain and depends only on N 
and D 0 but not on q and X. The value D eff allows us to 
obtain the number of particles per chain N. In the op- 
posite limit (i.e. l/a< 1), we are sensitive to motions on 
the size of individual particles. The main contributions 
to the measured diffusion coefficient come from inter- 
nal motions of the chains (i.e. particles’ fluctuations 
inside the chain). However, the reason for a linear de- 
pendence of D t gOX\ q is not clear. 

The fact that the initial value of D ej /D 0 (here equal 
to 0. 1 6 when q approaches zero) docs not depend on X 
indicates that the straight chain configuration is well 
maintained even for the lowest X values we used. This 
value can be used to find the number of particles per 
chain from the equation (4). This gives us N=19±l as 
the mean number of particles in a chain (after 6 hours 
at A.=406). 

Figure 6 shows the dependence of D e g on X based 
on the data in Fig.5. Since D e /D 0 ~k(qa), we plotted, in 
a Log-Log plot, the slope k versus X in Fie.6. Here, a 
fit of power-law is obtained D e /D 0 x X giving the 
exponent 5 =- 0.47 ± 0.05. 



We can approximate D ,.,/D 0 ~ 1/V X . We can see 
here that when the magnetic field increases (so X in- 
creases) this reduces the internal motions. Therefore, 
based on the results of Fig.5 and 6, may be 

written as: 


D eff /D 0 = — g(N)qa + f(N) (6) 


VI 


This effective diffusion coefficient given by Eq.6 is 
consisted of two parts. 

The first term is the contribution due to internal mo- 
tions of the chain (particle fluctuations) while the sec- 
ond part is due to the whole chain motion (diffusion of 
the center of mass of the chain). 

In Eq.6, g(N) accounts for the chain length depend- 
ence. g(N) and the constant A is determined from ex- 
periments described below. 

We created two chain lengths to study the effect of 
chain length on the effective diffusion coefficient. To 
obtain a longer chain size than 19 particles per chain, 
we hold the magnetic field at X=406 for 3 more hours 
(this increases the experimental time up to more than 9 
hours). We then reduced the field down to X=46 and 
repeated the q-dependence measurements. The result is 
shown in Fig.7 below where for the same value of X 
(X=46) two experiments with two different chain sizes 
are compared here: 


Figure 7: Effect of different chain sizes on the meas- 



From the data shown in Fig.7, we notice that when 
N varies, not only the initial value is different but also 
the slope changes. If we also assume that D t /D 0 °c 
frf ,e is found to (-0.7±0.1). 

To summarize the results from figures 5,6 and 7, 
we found that the effective diffusion coefficient fol- 
lows the simple relation: 

D - ,D ° = i^7r o+/(A ' , <7) 

Here f(N) is given by equation 4. The constant A 
was found equal to 3.0+0.2. It was determined experi- 
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mentally by Finding the exact value of D e gfDo from 
Fig. 5 for a given qa/27t. 

It is worth noting here that our experiments do not 
cover the whole range of probing window (l-2n fq). 
We were only sensitive to the most interesting range 
where a q dependency of D e % can be observed. How- 
ever, if we expand the q range toward both ends, dif- 
ferent behavior of D eff may be observed. Based on dy- 
namics of polymers [14,15], we expect a qualitative 
behavior as shown in the Fig. 8. Here, three dynamical 
regimes exist. 

For small qa values, D e /D p corresponds to the dif- 
fusion coefficient of the entire chain. In the opposite 
limit {qa>2n ) it corresponds to the diffusion coeffi- 
cient of the monomer inside a polymer chain D p . 



Figure 8: Universal behavior of the effective diffusion 
coefficient versus q in polymer physics. 


In this work, we just reached the boundary of the 
two limits. If we extend qa a little large, we may see 
the saturation behavior. From equations (7, 8 and 9) 
we can see that the higher X and N arc, the lower the 
internal vibration frequency and the slower the diffu- 
sion motion. Increasing the magnetic interaction leads 
to a more rigid chain, which is reasonable. Decreasing 
X reduces the dipolar binding between particles in a 
chain. Thus, fluctuations of higher frequency are ob- 
served. However, X can not be reduced below the 
value of 1. In this case, only the motion of isolated 
particles will be probed (if no irreversible aggregation 
occurred) since chains break [16]. The N dependence 
indicates that when the number of particles per chain 
increases, internal frequencies of the chain reduce. 
This behavior may be explained by the fact that more 
vibrational modes are present for a longer chain than 
that of a shorter chain. If the slowest modes (or the 
longest wavelength) contribute mainly to the probed 
fluctuations, the characteristic frequency will decrease 
with chain length. 

For N>5,f(N)~ (3/2)N 0 75 . Therefore, we may re- 
place AT 07 in Eq.6-9 by f(N). Thus, equation (7) be- 
comes: 


D eJf =(—p=qa + 1)/ (N)D 0 


'VI 


( 10 ) 


We can write Eq.10 by separating the two parts. The 
internal part with an apparent internal diffusion coeffi- 
cient D inh and the global part with the diffusion coeffi- 
cient of the whole chain D chain . 


In our case D p is the diffusion coefficient of a sin- 
gle particle in the chain configuration. Notice that D p 
is different from the case of isolated particles D 0 . In 
both of these regions, the behavior of D e f/D p does not 
depend on qa/2n . If a broader range of qa/2n was 
accessible we expect to see also these two others re- 
gions where D e j/D 0 becomes constant. 


In our experiment, the characteristic frequency [15] 
probed by DLS is defined as Q (q)=D eff .q 2 . 


«(?) = 


N 0 7 VI 


qa+f(N ) 


q 2 D o 


( 8 ) 


The two constant regions can be obtained by setting 
qa=0 and qa=2n . When qa=0, D eff =D cham or, in 
terms of the minimum frequency, Q min~~Q F) chain- 
When qa/2n approaches 1, the maximum frequency 
is obtained by replacing qa by 27t in Eq.7 for the ef- 
fective diffusion coefficient. The result, valid for 
q>(2n / a ), is: 


( - — + f (N))q 2 D 0 (9) 

N 01 VI 


D e ff - £*int + & chain 

With 

Djm = ~J=’D cllain and D -/(ff)O. 

VI 

Equations (8) and (9) can also be expressed in a 
closed form: 


n(q) = (-^q 3 a + q 2 )f(N)D 0 00 


VI 

And, for q>(2n /a): 

Q lm% (q) = {^+\)f(N)q 1 D 0 02) 

VI 


The apparent diffusion coefficient of particles in a 
chain configuration D p is given in terms of the maxi- 
mum frequency (term “apparent” is used since a parti- 
cle does not have a real diffusion behavior inside a 
chain because of its bounding motion): 

D p =- =(-£ + 1 )D cham (13) 

q yjA 

From equation (13) we see that D p reduces to D chain 
when A»(47t) 2 =160 and that the particle radius is 
completely contained in the single particle diffusion 
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coefficient D 0 . Because of the square root dependence, 
even in the case k=10 4 , internal fluctuations still con- 
tribute f 12% of that of the whole chain. To reach 
D in /D chain = 1 %, X=\0 6 is necessary. In practice, such a 
large A. value is difficult to obtain with submicronic 
particles, due to saturation of magnetization. In our 
case, if no saturation effect is taken account, the mag- 
netic field would have the value of 1 Tesla to get a 
completely rigid chain. In our case, for the maximum 
magnetic field used, A*=406 and D m /D chain - 62%. 

Physically, this means that in almost all experi- 
ments realized with these systems composed of submi- 
cronic super-paramagnetic panicles, chain fluctuations 
can not be neglected. Furthermore, the only way to 
obtain rigid chains is to increase the particle size well 
above the micrometer. Since the magnetic interactions 
is proportional to the cube of the particle radius (cf. 
Eq.l), a particle diameter about 3 micrometers need to 
be used with a magnetic field of 300 Gauss to reach 
?i=10 6 (assuming the same magnetic susceptibility). 
The biggest problem with such particles will be their 
sedimentation and stability of ferro fluid droplets. Fur- 
thermore, using these particles we can test if rigid 
chains coalesce into bigger aggregate as theoretical 
prediction suggested [17]. This is our next project. 

Conclusion 

In conclusion, we find that DLS is a powerful tech- 
nique in probing the dynamics of chains. We made 
quantitative measurements concerning the variation of 
our parameters; the coupling parameter X, the mean 
number of particles per chain N and the particle radius 
a. For instance, we showed that the frequency of the 
fluctuations depend on the magnetic interaction 
through I/V X (the higher the field, the more rigid 
the chain). The dependence of the number of particle 
per chain also leads to slow down the characteristic 
frequency when the chain length increases. 

We found separated motions when the scattering 
wave vector varies. The effective diffusion coefficient 
comprises two parts. The first one comes from the dif- 
fusion of the center of mass of the chain (drift motion), 
and the second one comes from the internal fluctua- 
tions of particles that compose the chain. These inter- 
nal fluctuations can not be neglected even for the 
highest magnetic field used in our experiments 
(k= 406, D in /D chain = 62%). This work is a basic study of 
the chain dynamics. It is the first step leading to a bet- 
ter understanding the mechanism of chain-chain lateral 
aggregation in the field-induced structural formation of 
dipolar fluids. 

Acknowledgements 

We gratefully acknowledge M. Hagenbuechle for 
the experimental setup and earlier work [9,18] on 
which this work is based and the NASA for the grant 
NAG 3-1830 supporting this work. 


References 

[1.] R.E. Rosensweig, Ferrohydrodynamics, Cam- 
bridge University Press, 1985. 

[2.] P.G. de Gennes & P.A Pincus. (1970), Phys Kon- 
dens Mat., 11, 189. 

[3.] G.A. Flores, J. Liu, M. Mohebi & N. Jamasbi, 
“Field induced columnar and bent-wall-like patterns in 
a ferrofluid emulsion”, proceeding of the 6th interna- 
tional conference on E.R.F . & M.R.S &A.T. Japan 
(1997) To be published. 

[4.] J. Liu et al, “Field-Induced Structures in Ferrofluid 
Emulsions”, Phys. Rev. Letters , 74 (14), p. 2828 
(1995). 

[5.] S. Cutillas and G. Bossis, “A comparison between 
flow-induced structures in electro- and magneto- 
rheological fluids”, Europhys. Lett ., 40 (4), 465 
(1997). 

[6.] E. Lemaire et al, “Influence of the particle size on 
the rheology of magnetic suspensions”, J of Rheology , 
39(5), (1995) 

[7.] S. Cutillas, G. Bossis and A. Cebers, ‘Flow- 
Induced transition from cylindrical to layered patterns 
in MR fluids” Phys Rev E, 57 (1), 804, (1998). 

[8.] B.J. Berne & R. Pecora, Dynamic Light Scatter- 
ing, Wiley and Sons, New York, (1966). 

[9.] M. Hagenbuechle and J. Liu, “Study of chain for- 
mation and chain dynamics in a dilute magneto- 
rheological fluid”, Applied Optics , 36, p. 7664, (1997). 
[10.] K. Zhan, R. Lcnke and G. Maret, J. Phys. II 
France , 4, p. 555, (1994) 

[11.] H. Yamakawa, Transport properties of polymer 
chains in dilute solution: hydrodynamic interactions, J. 
Chem. Phys., 53, p.436, (1970). 

[12.] G. K. Batchelor, “Slendcr-body theory for parti- 
cles with arbitrary cross-section in Stokes flow”, J. 
Fluid. Mech., 44, p.419, (1970). 

[13.] M. Fermigicr & A.P. Gast, “Structure Evolution 
in a Paramagnetic Latex Suspension”, J. Coll. Int. Sci ., 
154, 522,(1992). 

[14.] R. Pecora, Dynamic Light Scattering, Plenum 
Press, New York (1985). 

[15.] P. G. de Gennes, Scaling concepts in polymer 
physics, Cornell University Press (1979) 

[16.] S. Cutillas, G. Bossis, E. Lemaire, A. Meunier 
and A. Cebers, Experimental and theoretical study of 
the field-induced phase separation in electro- and 
magnetorheological suspensions”, proceeding of the 
6th international conference on E.R.F. & M.R.S ScA.T. 
Japan (1997) To be published. 

[17.] T. C. Halsey and W. Toor, J. of Stat. Phys., 61, 1257, 
(1990) 

[18.] M. Hagenbuechle and J. Liu, proceeding of the 
6th international conference on E.R.F. & M.R.S &A.T. 
Japan (1997) To be published. 


105 


